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I Using Weierstrassian elliptic functions the exact geodesies in the 

£N| ■ Schwarzschild metric are expressed in a simple and most transparent 

form. The results are useful for analytical and numerical applications. 
For example we calculate the perihelion precession and the light deflec- 
tion in the post-Einsteinian approximation. The bounded orbits are 
computed in the post-Newtonian order. As a topical application we cal- 
culate the gravitational red shift for a star moving in the Schwarzschild 
field. 
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1 Introduction 



Schwarzschild geodesies are elliptic functions, therefore, they should be writ- 
ten as such. For this purpose the Weierstrassian elliptic functions are most 
useful because they lead to simple expressions. The reason for this is that 
the solution of quartic or cubic equations can be avoided in this way. 

In a recent paper [1] an analytic solution for the geodesic in the weak- 
field approximation was given. As pointed out in that paper the progress 
in the astronomical observations call for better analytical methods. In this 
respect it is desirable to have the exact geodesies in a form most suited for 
applications. For the orbits in polar coordinates (next section) this goal can 
be achieved by using Weierstrass' 'P-function for which many analytical and 
numerical methods are known [2]. Considering the motion in time (section 
4) the related (- and u-functions of Weierstrass appear. 

Jacobian elliptic functions have been used by Darwin [8] for the form of 
the orbits. After some transformation our result (2.13) agrees with his. But 
in his second paper he abandons the elliptic functions because they were 
"not so well adapted to a study of the time in those orbits" . Obviously the 
Weierstrass functions are better suited for the problem. Indeed, expressing 
them by Theta functions one gets the natural expansion of the geodesies 
in powers of the Schwarzschild radius, this expansion involves elementary 
functions only. The Weierstrass functions have also been used by Hagihara 
[11] 1 . But he has chosen the variables and constants of integration in a 
manner which leads to less explicit results. So it is difficult to derive the 
post-Newtonian corrections to the geodesies given here from his formulas. 
As a topical application we finally calculate the red shift for a star moving 
in the Schwarzschild field. The geodesies are also needed for the study of 
modifications of general relativity ([10], section 5.13). 

2 The orbits in polar coordinates r = r((f) 

We take the coordinates x° = ct, x 1 = r, x 2 = i?, x 3 = ip and write the 
Schwarzschild metric in the form 

ds 2 = T —^c 2 dt 2 — dr 2 - r 2 (d$ 2 + sin 2 -ddip 2 ) (2.1) 

r r — r s 

1 1 am indebted to C. Lammerzahl and P. Fiziev for bringing this reference to my 
attention 
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where r s = J* is the Schwarzschild radius. We shall assume c = 1 in the 
following. The geodesic equation 

d 2 x a ^„ dx" dx 1 

with the Christoffel leads to the following three differential equations 

d 2 t fdt dr 
ds 2 ds ds 

Here we have used the standard representation 

1 - - = e v (2.6) 

r 

and have chosen # = ir/2 as the plain of motion. The Christoffel symbols 
can be taken from the Appendix of [3]. 
Multiplying (2.3) by expz/ we find 



ds ( ds ) 



ds v ds' 

so that 



,.dt 

e — = const. = E 
ds 

Ee~ v . (2.7) 



ds 

Next multiplying (2.5) by r 2 we get 

2 dt f . r 

r — = const. = L, 

ds 

hence 

dip L 
ds r 2 

For the constants of integration we use the notation of Chandrasekhar [4] . 



(2.8) 
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Finally, substituting (2.7) and (2.8) into (2.5) and multiplying by 2 exp(— v) 
xdr/ds we obtain 



d_ 

ds 



-»(±) 2 -E*e-» + ^}=0. (2.9) 



Consequently, the square bracket is equal to another constant = b. Then 
the resulting differential equation can be written as 

The constant b can be arbitrarily adjusted by rescaling the affine parameter 
s. Below we shall take b = —to 2 where m is the rest mass of the test particle. 
This will enable us to include null geodesies (light rays) with m? = 0. Each 
geodesic is characterized by two constants of the motion: energy E and 
angular momentum L. 

Taking the square root of (2.10) and dividing by (2.8) we get 



dr IE 2 — to, 2 . m? 



r 2 + r s r=JJ(r). (2.11) 



Now ip = <p(r) can be written as an elliptic integral. However, it is better to 
consider the inverse r = r((f) in terms of elliptic function by using a formula 
of Weierstrass ([5], p. 452). Let the quartic f(r) be written as 

f(r) = a r 4 + 4air 3 + 6a 2 r 2 + 4a 3 r + a 4 , (2.12) 

and let T\ be a zero f{r\) = 0, then a solution of (2.11) is given by 

r = n + AV(^g 2 f g { 3 )-f"(n)/6- (2 ' 13) 
Here V((p;g2,g3) is Weierstrass' P-function with invariants 

<?2 = oo«4 — 4aia3 + 3a 2 . (2-14) 

#3 = aoa2a4 + 2aia2a3 — a\ — aoa| — a\a^. (2.15) 

In our case we have 04 = 0. For the convenience of the reader we reproduce 
the short proof in the Appendix. 

The result (2.13) is not yet the solution of our problem because it con- 
tains too many constants: the invariants g 2 ,g3 and the derivatives of / 
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depend on E, L, but in addition the zero r\ appears. Of course one could 
calculate r\ as a function of E, L by solving the quartic equation f(r) = 0, 
but this gives complicated expressions. It is much better to use r\ and a 
second zero r 2 as constants of integration instead of E, L. This is even de- 
sirable from the astronomers point of view because the zeros of derivative 
(2.11) are turning points of the geodesic, for example in case of a bounded 
orbit they can be identified with the perihelion and aphelion of the orbit. 
In order to express E, L by r\ , r 2 we write our quartic in the form 

f(r) = a r(r - n)(r - r 2 )(r - r 3 ) (2.16) 

and compare the coefficients of r 3 ,r 2 ,r with (2.12). This leads to 

2 

171 

4ai = -a (n + r 2 + r 3 ) = jjr s 

602 = ao(rir 2 + rir 3 + r 2 r 3 ) = -1 (2-17) 
4a 3 = -a nr 2 r 3 = r s . 



Since 



we can solve for 



E 2 - m 2 , 
a = L2 (2.18) 

%r s = r \ +r2+ ; 3 , (2.19) 

E 2 -m 2 = ^ . (2.20) 

r\ +r 2 + r 3 

In addition we obtain the third zero 

r 3 =r s ^ . (2.21) 

nr 2 - n r s - r 2 r s 

The relations (2.19-21) allow to express everything in terms of r±,r 2 . For 
the invariants we find 

= J_ _ Z^ r 2 = _L _ r s n+r 2 + r 3 22 
52 12 4L 2rs 12 4 nr 2 + nr 3 + r 2 r 3 1 ' J 

1 m 2 2 oo 2 
93 = ¥ ~ 48l2 rs " 16 rs 
1 1 (ri + r 2 + r 3 )r s | 1 r 2 
6 3 48 r\r 2 + rir 3 + r 2 r 3 16 r\r 2 + rir 3 + r 2 r 3 ' 
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Here r 3 has to be substituted by (2.21). For the derivatives f'(ri), f"{r{) 
which appear in our solution (2.13) we obtain 

f( ri ) = - ri(ri ~ r2)(ri ~ r3) (2.24) 

r x r 2 + rir 3 + r 2 r 3 

fi/f ri \ = 2 ^ l - T2)ri + - r ^ n + - r2 ^ n - r ^ (2 25) 

nr 2 + rir 3 + r 2 r 3 

With these substitutions the result (2.13) gives all possible geodesies in the 
form r = r(tp; r±, r 2 )- This will be discussed in the next section. 

As a first check of the solution (2.13) we consider the Newtonian limit. 
Let the two zeros r±, r 2 be real and large compared to the Schwarzschild 
radius r s in absolute value. Then neglecting 0(r s ) in the invariants (2.22-23) 
the "P-function becomes elementary ([2],p.652, equation 18.12.27): 

Pte i 6 ~ 3) = ~h + 4^72 = -h + 2(1 -L v y (2 ' 26) 

The leading order in the derivatives of / is given by 

fin) = -— (n - r 2 ) 

f'(n) = -2(3^ -2). (2.27) 
\ r 2 > 

It is convenient to introduce the eccentricity z by 

r\ 1 — e 
r 2 1 + e ' 

Using all this in (2.13) we find the wellknown conic 



(2.28) 



r= (1 + g)n . (2.29) 

1 + £ COS if 

Assuming both zeros ri, r 2 positive and r\ < r 2 we have e < 1 and the orbit 
is an ellipse with perihelion n and aphelion r 2 . In the hyperbolic case e > 1 
we see from (2.28) that if r\ is positive r 2 must be negative. Then there is 
only one physical turning point r± which is the point of closest approach. 
The latter always corresponds to (p = 0. The relativistic corrections to (2.29) 
are calculated in the following section. 
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3 Discussion of the solution 



The solution r = r(<p) (2.13) is an elliptic function of (p which implies that it 
is doubly-periodic ([2], p. 629 or any book on elliptic functions). The values 
of the two half-periods depend on the three roots of the fundamental 
cubic equation 

4e 3 - g 2 e - g 3 = 0. (3.1) 

Again it is not necessary to solve this equation because the solutions ej,j = 
1, 2, 3 can be easily obtained from the roots 0, r±, r 2 , r% of our quartic f(r) = 
0. To see this we transform /(r) to Weierstrass' normal form as follows. 
First we set r = l/x so that from (2.16) we get 

f(r) = -^(4a3X 3 + 6a 2 x 2 + Aa\x + ao) 



Next we remove the quadratic term by introducing 

i =x = i(e-|). (3.2) 
r as v 2 / 

This gives the normal form of Weierstrass 

/(r) = 7 1 , 4 (4e 3 -g 2 e- g3 ), (3.3) 
(e- |a 2 ) 4 

with the above invariants (2.14-15). That means roots of /(r) are simply 
related to roots of (3.1) by the transformation 

e , = ^ + ^ = ^-l. (3.4) 
3 rj 2 Arj 12 V ; 

The cubic equation (3.1) with real coefficients has either three real roots 
or one real and two complex conjugated roots. The first case occurs if the 
discriminant 

A = gl - TJg\ (3.5) 

is positive, in the second case A is negative. In terms of the roots A is given 
by ([2], p. 629, equation 18.1.8) 

A = 16(ei - e 2 ) 2 (e 2 - e 3 ) 2 (e 3 - exf . (3.6) 

The physically interesting orbits correspond to the first case of real roots. 
If we have two complex conjugated zeros r 2 = r\ then (2.28) implies that 
the eccentricity e is imaginary. Such orbits have been discussed by Chan- 
drasekhar ([4], p. 111). Now we discuss the various cases. 
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3.1 Bound orbits 



In this case we have two positive turning points r 2 > r\ > 0, consequently 
there are three real roots e\ > > e 2 > e 3 given by 



ei 



r s 1 1 r s n + r 2 
4r 3 12 ~~ 6 4 nr 2 

1 r. 



e 2 



e3 



(3.7) 
(3.8) 



12 ' An 7 12 + 4r^' 

Our convention is chosen in agreement with [2]. The real half-period u; of 
the P-function is given by ([2], p. 549, equation 18.9.8) 



dt 



K(k 2 



^4t 3 -g 2 t-g 3 y/d ~ e 3 



(3.9) 



where K(k 2 ) is the complete elliptic integral of the first kind with parameter 



k 2 = ^^=r 



r 2 - n 



e 1 -e 3 nr 2 - r s (2r 1 + r 2 ) 
r s 2e / 3-ey 



(3.10) 



3-e\-i 

ri i+£v i_rs rn 

As a first application let us give the post-Einsteinian correction to the 
orbital precession. If k 2 (3.10) is small we can use the expansion ([2], p. 591, 
equation 17.3.11) 



7T 



1\2, 



1 • 3\2 



2-4 



fc 4 + ... 



(3.11) 



From the roots e\, e% we find 



= 2 



1 + 



r s 2n + r 2 3 2 / 2r i + r 2 



This finally leads to the half-period 



H — r 

2 rir 2 8 s 



nr 2 



+0(r 



r 3r s n +r 2 

CJ = 7T|1 H ; h 



3rf 



4 nr 2 ' 8(rir 2 ) 2 
3r, 1 



+ 



(2n + r 2 )(2r 2 + n) + -(r 2 - n) 2 ] } 
3 r 2 18 + e 2 



nl + e 18r 2 (l + e) 2 



+ 



0(r s 3 )} 



(3.12) 
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The perihelion precession is given by Aip = 2(lj — ir). Then the order r s 
in (3.12) is Einstein's result and the 0(r 2 ) gives the correction to it. The 
accurate computation of the half-period is necessary to control the orbit in 
the large. 

To compute the relativistic corrections for r((p) from (2.13) we express 
the "P-function by Theta functions ([5], p. 464) 

$i(z, q) = 2(/ 1 / 4 (sinz — q 2 sin3z + q 6 sin5z — . . .) 

$2{z-,q) = 2(7 1 / 4 (cosz + g 2 cos3z + q 6 cos5z + . . .) 
#3(2,9) = 1 + 2g(cos2z + <7 3 cos4z + q 8 cos6z + . . .) 
•&i(z,q) = 1 — 2g(cos2z — q 3 cos Az + q 8 cos6z —...). 
Here q is the so-called Nome ([2], eq. 17.3.21) 



(3.13) 



k 2 

q= Te + 



(3.14) 



These series are rapidly converging since k 2 is small (3.10), they give the 
natural expansion in powers of the Schwarzschild radius r s . Now the V- 
function is given in terms of Theta functions by ([2], eq. 18.10.5) 



e 2 + AiO 2{ M0)M(j)) ) 



(3.15) 



where 



Using 



1 f s 

"12 4r7 



7T 



ALU 2 



sin 



1 + 4g(cos 2 <f>-l) + 0{q 2 ) 



this leads to 

f(n) 



V{<p)- 



/'(ri) = 2n( T 

f"(n) = 



r, 3e + e 



+ e n(l + £f 



1 - 5e 



+ 6 



. r, 1 — Ae — e 2 



l + £ n (l + e) 2 



24 



l + e - 2esin 2 < 
4(1 + £) sin 2 4> 



(3.16) 



(3.17) 



1 + 



r\ 1 + £ cos 2(f) 



(-3-|(l-cos0)+ 
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+2e i±£ s in 2 ^)]+0(^) . (3.18) 

Substituting this into (2.13) gives the desired orbit to 0(r s ) 
r(ip) l + e r s 2sin 2 



r s lsm"<p r 1 £, . 

H £ ~-> 777 \ --, 777 3 + 7:(1 _ COS0)- 

n l + ecos2(/>ll + ecos2(/>L 2 V ^ ; 



r\ l + e cos 2^ n l + e cos 20 ^ 1 + e cos 20 

_ 2e ^±£ s in 2 0l-^}. (3.19) 
l + e J 1 + eJ v ' 

It is important to insert the period uj in (3.16) according to (3.12) in order 

to describe the perihelion precession correctly. 

If the two roots r\ = r 2 coincide, it follows from (2.24) that f'(ri) = 0. 

According to (2.13) we then have circular motion r = r\. If all three zeros 

coincide r\ = r2 = then (2.21)gives r% = 3r s which is the innermost 

circular orbit. 

3.2 Unbound orbits 

In this case there is only one physical point, the point of closest approach r\. 
The other root r2 is negative, therefore, it is better to use the eccentricity e 
(2.27) as the second basic quantity. With r% given by (2.21) we then have 

ri > r 3 > > r 2 = Y~~ ri ' ( 3 - 20 ) 

because e > 1. The periodicity of (2.13) in (p is now realized by a jump to an 
unphysical branch with r < 0. In reality a comet moves on one branch only, 
but it is a tricky problem to decide on which one. This is due to the fact 
that the period differs a little from 2ir as in the bounded case. Consequently, 
neighboring physical branches r > are a little rotated against each other 
and the distinction between them is not easy. The quantity of physical 
interest is the direction of the asymptote. It follows from the original 
equation (2.11) by integrating the inverse over r from n to oo 

oo 

dr 



/dv 
jm- (3 - 21) 



n 



This is an elliptic integral which can be transformed to Legendre's normal 
form 

a f d<& , 
^ = -^=\ , (3.22) 



^1 - k 2 sin 2 4> 
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by the transformation ([6], vol.11, p. 308) 



S in^= r -^l r -^±, sin 2 <t> 2 = ^^. (3.23) 

T\ — r 2 T — T3 T\ — r 2 

The parameter k 2 in (3.22) is given by 

fe2 = rsn-r,^ 
r\ r 3 - r 2 

and 

H = 2 ■ (3.25) 

Vri(r 3 - r 2 ) 

The integral (3.22) is an incomplete elliptic integral of the first kind 

^ OD = ^F($ 2 ,k 2 ) (3.26) 



which has the expansion ([6], vol.11, p. 313) 

k 2 1 

F($ 2 , /c 2 ) = $ 2 + — ($ 2 - - sin2$ 2 ) + 0(k 4 ). (3.27) 
For small r s /Vi we find 

k 2 = ^ r ^ + o( r ^) 2 

e + 1 r\ Vri/ 

l + eri 
sm $ 2 = ^r— 1 ) • 



This gives 



and 



1 £ ~ 1 r 3 
cos 2$2 = 

e e r\ 



^ 00 = 2$ 2 + ^(^-$ 2 --^-sin2$ 2 ). (3.28) 
ri \e + 1 2e + 2 / 



?£/ 3 _ 

n 

It is convenient to calculate 



The leading order is the Newtonian asymptote of the hyperbola. 
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3.3 Null geodesies 

For m? = there is only one constant of integration in the quartic (2.11) 

4 

which is the so-called impact parameter 

d=|. (3.30) 

Now it is necessary to calculate the roots of f(r) = 0. This is easily done 
by means of a power series expansion 



We find 



r = c d + c\r s + c 2 r s + . . . 
5 3 



n = <l~-^ 2 + 0(* 



r 2 = -d(l+ 5 --h 2 + 0{5')) (3.31) 



2 8 
r3 = d(<5 + ( 5 3 ), 



where 



5 = 1 ( ' 3 ' 32 ' 1 

and we have ordered the zeros in the same way as in (3.20). Then as in the 
last subsection we can calculate the direction of the asymptote (3.21) which 
now is equal to 

< Poo = ndF^k 2 ) (3.33) 

with $i given by (3.23) 

sin 2 $ 1 = _ + _$ + 0(5 3 ) (3.34) 

and k 2 by (3.24) 

k 2 = 25{l-5 + ^5 2 ) (3.35) 

and fi by (3.25) 

(3.36) 
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We want to calculate the light deflection in the post-Einsteinian approx- 
imation. Using again the expansion (3.27) we have 

^oo = 2(l - 5 - + j^ 2 ) [<&! + ^(1 - <5)($! - \ SU12*!) 



From (3.34) we obtain 



Then up to 0{5 2 ) we find 



7T . /3 3 \ r o 



The deflection angle is given by 



(3.37) 



(3.38) 



Instead of the impact parameter d in 5 (3.32) we would like to use the 
distance of closest approach n (2.24) in the form 



Si 



n 



The two are related by 



which leads to 



S = 6l --S 2 1 + 0(6l) 



A^ = 25 1 + (i + ^)<5 2 . 



The first term 25\ is Einstein's result. 



(3.39) 



(3.40) 



4 The motion in time t = t((p) 

By dividing (2.7) by (2.10) we find 

E r 3 
~~L(r-r 3 )y/f(ry 



dt t 

T = E 

ar r — r, 



-1/2 



(4.1) 
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We choose t = at the point of closest approach r = r\ and get 



E 

T 



dx 



n 



v/7C 



X 



x 2 + r s x + r 2 s + 



:)■ 



(4.2) 



This is a sum of elliptic integrals of first, second and third kind. Calculating 
these gives the coordinate time as a function of r. However, we want t as 
a function of ip and, therefore, use the substitution (6.7) of the Appendix 
again 

/'(n) 



x = ri + 



4P(<p) - /"( ri )/6 

dx 

dtp, 



(4.3) 



v/7R 

where the last relation follows from (2.11). 

The last integral 0(rf) in (4.2) is a small correction and we neglect it at 
the moment. Then integrals of the following form remain to be calculated 



J, 



du 



(V(u) - V(v)Y 



(4.4) 



where we have set 



V{v) 



fin) 



24 



Such integrals are known ([7], vol.4, p. 109-110) 



V'(v) 



2C(w)(^ + log 



a(v — ip) 



<j(v + ip) 



M<p) 



V' 2 {v) 



(4.5) 



(4.6) 



C(ip + v) + ((ip -v) + 2V(v)p + V"(v) J^tp) . (4.7) 



These results are easily verified by differentiating and using addition for- 
mulas. Of course Jo(<p) is just the polar angle ip. Then (4.2) leads to the 
desired result for t(ip): 



E 

= j{ w + nJiM + T 2 J 2 {^)}+0{rl). 



where 



tq = t\ + r s ri + r 2 s 



(4.8) 
(4.9) 
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n = (y + j)/'(n) (4-10) 
fiii 



* - ££> (111) 

Again we evaluate this for bounded orbits in the post-Newtonian approxi- 
mation by means of the expansion in Theta functions. 

The quantity v in (4.8-11) is given as the zero of (3.18). Introducing 

V = — v (4.12) 

we find 

cos2 y = -I-!£^±i-^- / 3. (4.13) 

e n 2s v ' 

Since e < 1, V is complex: 



2y = 7r + ilog(/3 + ^//3 2 - l)= Tr + 2ib (4.14) 
Using ([2], eq.18.10.6) 

vr 3 MV)Mv)MvWHo) 



r'(v) = 



4w 3 i? 2 (0)t? 3 (0)i?4(0)<?f(y) 

7T 3 COSU , 

-4^^7 + ° ((?) 



we obtain 



' W 2 W 3 (/3 + 1) 2 + [ s) 



= ^(TTiF( 1 + n i-e 2 )■ (415) 

Similarly we calculate ((v) from ([2], eq.18.10.7) 
where 



and cr(z) from ([2], eq. 18. 10.8) 
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This implies 

a(v — (p) n , s'm(V — <&) ^. 9 . 

log ^ = -2-V + log . + 0(q 2 ) 
a{v + (p) uj sm(l/ + <P) 

= -2-V + log ;., ■ - = -2^v<p + 2ia(<p), (4.19) 

uj cos(z6 + $) uj 

where a is given by 

a(ip) = arctan(tan<I>tanh&). (4.20) 
Then J\ (tp) (4.6) is equal to 

To expand this in the post-Newtonian order we first calculate a(<p) from 



1-e x r s l + 3e /1-e 
tana = \ tan5> H e K \ tan5>. 



Introducing 



V> = 2arctan(y / ^-^tan$) (4.21) 

we get 

ib r s e 1 + 3e . , 

a = — + — ^ sm ^- 

2 n 4 1 — e l 

As before in (3.19) we do not expand in (4.21). However, if one does so 
one finds a contribution 0(r s jr\) 

ib = 2arctan(j— ^ tan(w/2) N ) -- — W — (4.22) 

r y\]l + e ') 2ri\/ l + el + ecosv? v y 

where the first term, say ipN, is the parameter which appears in Newtonian 
mechanics (Kepler's equation, see below (4.29)). Now the expansion of J\{ip) 
is given by 

J 1 {<p) = J?{<p) + ^Jl(<p) 



J°( v ) = 2^(J±±±-<p) (4.24) 
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-<p(7 + 3e)+e 



l + 3e /l + e 



' 1 + e 5 - 4e - 3e 



2 n 



, l-e l-e 
For e — )■ we have the simple finite limit 

Ji(p)| e =o = — (l + 3-) (2$ -sin 2$) 



7T V 



which also follows directly from the definition (4.4). 
To expand 



M<p) = 



V' 2 (v) 



,f] „/ 7T sin 2$ ... . . . 



we need 



J 2 V) + -ij 2 V) 



7?"(u) = 6V 2 (v) - § 



e 2e - 1 



r, 11 - 18e - 5e" 



2(l + e) 2 V 2n (2e- l)(l + e 



■)■ 



This finally leads to 



= - 



4 (l + e) 



e 1 



95 ■ 



(l + e) sin 2$ 2e) - 1 



l + ecos2$ 2(1 + e; 



2 (l + e) 



e 1 — e 



15 - 12e - lie 2 



sin2<I> 



■( 



1 + 4e + 3e z 



1 + e cos 2$ V 1 + e cos 2<3? 



+ 



l-e 2 

17 - 13e - 12e 2 
T^e 



+ 



o 6-15e + 8e 2 -e 3 x 1 - 2e 



(4.25) 



(l + e) 2 (l-e) "—l + e J" (4>26) 

Again we do not expand <I> (3.16) in order to keep the perihelion precession 
as precise as possible. The limit for e — > is equal to 

J 2 (^)U=o = 2^(l + 6^)(6$ - 4 sin 2$ + sin 2$ cos 2$). 
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In the final result (4.8) for the time 



t(ip) =t (¥>) + -ti(¥>) (4.27) 
n 



the pre-factor £7/L also gives a correction: 



L Yr s ri(l + e)V n(l + e)/ 
which follows from (2.19-21). In to the terms proportional to ip cancel 



/ 2r? 1+e/ ib esin2$ \ 

Approximating 2$ by ip this is in agreement with Kepler's equation 



The post-Newtonian corrections in (4.27) come from various places. To show 
this we write the result in the form 



l + e \r s (l+e) 



1 + e 



As in (3.19) the post-Newtonian correction vanishes for circular motion e = 
0. 

5 Gravitational red shift 

The study of Schwarzschild geodesies is relevant for the investigation of the 
recently discovered S-stars near the Galactic Center ( [1] and references given 
there) . These stars move in the strong gravitational field of the central black 
hole so that general relativistic effects are observable and the Schwarzschild 
metric g^ u is a fairly good description of the situation. The measurable 
quantity of interest is the red shift of spectral lines in the light emitted by 
the moving star. Therefore we finally consider this. 
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Let v\ be the frequency of a given atomic line from the star and uq the 
frequency of the same line observed in the rest frame of the galaxy. If dx^ / dt 
is the velocity of the star, the two frequencies are related by ([9] p. 83, equ. 
3.5.6) 

We assume that the observer at X is far away from the center such that 
the denominator can be approximated by 1. For a star moving in the plane 
■d = it/ 2 we have 

. .dx a dx 13 ,. „{dr\2 9/ <d(p\2 ,„ 

«WjT-'-'(l) - r W ■ < 5 ' 2 > 

From (2.10) and (2.7) we find 



and (2.8) gives 

dip _ Le^_ 
~dt ~ Er^' 

Substituting all this into (5.1) we see that L drops out and we end up with 
the simple result 

— = ^ = — (1-— M. (5.3 

vq E E\ r(<f)J v ' 

By (2.20) we can express E by the perihelion n and aphelion r 2 

?2 2(-i r s \ 



E 2 =m 2 (l 



r\ + r 2 + r 3 / 

where r% is the small correction (2.21). Then we finally get 

vi /-, r s \ ( r s \-i/2 



^ = (l-Jj-)(l r JL r 1/Z . (5.4) 

The lowest order 0(r s ) is equal to 

l^-iti) + W^ + ^ (5 ' 5) 

Since the last term is always smaller than the second one we indeed have 
red shift Ai > Aq. Of course, it is maximal at the perihelion where r = r\ 
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is minimal. The total observed red shift is obtained by multiplying (5.5) 
with the Doppler factor (1 + v r )~ 1 where v r is the component of the relative 
velocity along the direction from the observer to the star ([9], p. 30). 
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6 Appendix: Integration of the differential equa- 
tion 

We closely follow Whittaker and Watson ([5], p. 452). With the notation of 
the paper (2.11), let 

} dx , . 

<P = -fW=, (6-1) 

where r\ is any zero, f(r\) = 0. By Taylor's theorem, we have 

f{x) = AA 3 {x - n) + 6A 2 (x - ri) 2 + 4A 1 (x - n) 3 + A (x - n) 4 , 
where 

A = a , A x = a ri + a\ 
A 2 = a Q r\ + 2a\r\ + a 2 , ^3 = ao^i + 3air 2 + 3a 2 ri + 03. 
Introducing the new integration variable 

s = (x-n)- 1 , s 1 = {r-r 1 )-\ (6.2) 

we have 

00 

<p = j[AA 3 s 3 + 6A 2 s 2 + 4A lS + A ]- 1/2 ds. 

si 

To remove the second term in the cubic we set 

and we get 

00 

(p = J [4z 3 - (3A 2 2 - AAxA^z - (2A 1 A 2 A 3 -A%- A Al)]-^ 2 dx (6.4) 
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The coefficients of z and z° are just the invariants 52, 53 (2.14-15) of the 
original quartic. 

Now the inversion of the integral gives Weierstrass' P-function 

zi=P((p;92,93)- (6-5) 

From (6.2) and (6.3) we have 

r = n + (6 - 6) 

and hence 

r = n + . (6.7) 

1 + 4P(<p) - /"(n)/6 l ° 
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